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Abstract. We consider asymptotics of ratios of random characteristic 
polynomials associated with orthogonal polynomial ensembles. Under some 
natural conditions on the measure in the definition of the orthogonal poly- 
nomial ensemble we establish a universality limit for these ratios. 



1. Introduction 

1.1. Formulation of the problem. In this article we consider orthogonal 
polynomial ensembles of n particles on K. Such ensembles are described by a 
positive Borel measure /i with finite moments, and the associated distribution 
function for the particles {x\, . . . , x n } of the form 

1 n 
dF^ n {x) = — A n (x) 2 Y]dfi(xi). 

Zn i=l 

Here Z n is the normalization constant, 

Z n = J ... J A n (x) 2 Y[dfi(xi), 



8=1 



and 

A n (*)= n 

n>i>j>l 

is the Vandermonde determinant. For symmetric functions f(x) = /(xi, . . . , x n ) 



[Xi Xj / 



of the x[s, 



[■ ■ [ f(x)A n (x) 2 f[d^ Xl ) 

Zjnj J i=l 



denotes the average of / with respect to <iP Min (x). 
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Given a set of n points, {xi, . . . , 



}, let 



n 




(a) = Y\_( a ~ 



i=l 



considered polynomial in a 6 C. When the chosen to be the 

random particle locations for an orthogonal polynomial ensemble, this becomes 
a random polynomial which we denote simply by D n (a). By extension from 
the random matrix case (see below), D n (a) is known as the characteristic 
polynomial for the orthogonal polynomial ensemble. 

The goal of the present paper is to study the large n asymptotics of the 
averages 



In particular, we focus on the bulk of the support for measures /x, and aim 
to establish universality of the scaling limit of averages (II. ip under mild as- 
sumptions on jM. Our main result says roughly that for ji locally absolutely 
continuous with a bounded Radon-Nikodym derivative, universality for the re- 
producing kernel implies universality for (jl.ip . It should be emphasized that 
aside for the existence of moments, we make no global assumptions on \i. In 
particular, we do not assume that \x is globally absolutely continuous or that 
supp(/i) is compact. 

1.2. Motivation and remarks on related works. In the case when d\x{x) = 
e~ v ^dx one can interpret the particles in the definition of the orthogonal poly- 
nomial ensembles as eigenvalues of a random Hermitian matrix taken from a 
Unitary Ensemble of Random Matrix Theory (RMT). The most well studied 
case is that of V(x) = x 2 , called the Gaussian Unitary Ensemble. In this 
context random characteristic polynomials are indeed the characteristic poly- 
nomials of a random Hermitian matrix. 

Averages of characteristic polynomials of random matrices are basic objects 
of interest in RMT, and were first considered by Andreev and Simons [1], Brezin 
and Hikami [B] and Keating and Snaith [13]. In particular, such averages are 
used to make predictions about zeros of the Riemann-zeta function on the crit- 
ical line (see, for example, Keating and Snaith [13], Conrey, Farmer, Keating, 
Rubinstein and Snaith [8], the survey article by Keating and Snaith [2], and 
references therein). Averages (II. ip are related to certain important distribu- 
tion functions studied in physics of quantum chaotic systems. Two examples 
involve the curvature distribution of energy levels of a chaotic system and the 
statistics of the local Green functions, in particular, the joint distribution of 
local denity of states; see Andreev and Simons [T] and references therein. Many 
other uses are described, for example, in Brezin [5]. 



(1.1) 



< 



.Dn(ai) ■ ■■Dni.Ctk) 

D n (/3x) . . . D n (/3 k ) 
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For averages ( 11. ip a number of algebraic and asymptotic results is available 
in the literature. Papers by Baik, Deift, and Strahov [3], Fyodorov and Strahov 
[12] and Borodin and Strahov [7] give explicit determinantal representations 
for fll.ip . These representations can be used for the asymptotic analysis as 
n — > oo. In particular, the asymptotics of (11. ip were investigated in Strahov 
and Fyodorov [12] in the case when dfx(x) = e~ v ^dx, and V(x) is an even 
polynomial. Strahov and Fyodorov [12J deal with the asymptotic in the bulk of 
the spectrum. Vanlessen [23] shows that for certain class of unitary ensembles 
of Hermitian matrices, averages (II. ip have universal asymptotic behavior at the 
origin of the spectrum. The asymptotic analysis in [12], and in [23] is based on 
the reformulation of an orthogonal polynomial problem as a Riemann-Hilbert 
problem by Fokas, Its and Kitaev [IT] . The Riemann-Hilbert problem is then 
analyzed asymptotically using the noncommutative steepest-descent method 
introduced by Deift and Zhou, see Deift [H] and references therein. 

In recent years, it has become evident that orthogonal polynomial ensembles 
play a role in probabilistic models other than RMT as well (see, for example, 
the survey paper by Konig [H]). As the relevant measures in these models 
are not necessarily of the form e~ v ^dx, it is of interest to study the prob- 
lem of universality limits for basic quantities of interest for a broader class of 
measures. Lubinsky's universality theorems regarding bulk universality for the 
reproducing kernel (see Lubinsky [161 Theorem 1.1], and [18| Theorem 1.1], and 
also Findley [10] , Simon [19], Totik [22] and Avila-Last-Simon [2] for exten- 
sions of Lubinsky's results and methods) are important steps in this direction. 
Our goal in this paper is to establish the corresponding universality limits for 
averages (II. ip in the bulk of the support of \i. 

1.3. Description of the main result. For n = 0, 1, 2, . . . we introduce the 
orthonormal polynomials associated with «, 

p n (x) = 7 n X n + • • • , 7n > 0. 

The orthonormality conditions are 



The nth reproducing kernel (also known as the Christoffel-Darboux kernel) 
for u is 



k=0 

If dfji{t) = w(t)dt in a neighborhood of x, we also define the normalized kernel 
to be 




n-l 




K n (x,y) = w(x)K n (x,y). 
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The Christoffel-Darboux formula enables one to rewrite K n (x,y) as 

/, x v , x ln-lPn(x)Pn-l(y) -Pn-l(x)Pn(y) 

(1.2) K n (x,y) = . 

7n x-y 

The reproducing kernel plays a special role in the theory of orthogonal poly- 
nomial ensembles. This is because the correlation functions for the orthogonal 
polynomial ensemble can be expressed as determinants of a matrix whose en- 
tries are given by values of the reproducing kernel (see, for example, Deift 

We say K n has a universal limit at x, if for any a, b G C we have 

/, ox r Kn(yX + ' X + Knfc)) sin7r(a-&) 

1.3 hm — r = — = %6. 

n-xx> K n (x,x) 7i(a — b) 

We say K n has a uniform universal limit at x, if the limit in ( 11. 3p is uniform 
for a, 6 in compact subsets of C. 

Two approaches introduced recently by Lubinsky ([IE]- [18]) make it possible 
to establish universality for K n under relatively mild conditions on //. These 
approaches were further extended and generalized by Findley [10], Simon [19] . 
Totik [22] and Avila-Last-Simon [2]. A typical result is: 

Theorem 1.1. Let [i be a probability measure on R with compact support that 
is regular in the sense of Stahl and Totik |20j. Suppose x G supp(yu) has a 
neighborhood, J, such that fi is absolutely continuous in J: dfi(t) = w(t)dt for 
t G J. Assume further, that w is positive and continuous at x. Then uniformly 
for a, b in compact subsets of the complex plane, we have 

K (x + — - x + — - ) 

1.4 hm £rf — r y -^ = S(a,b). 

In other words, K n has a uniform universal limit at x. 
Remarks. 

(a) Theorem 11.11 implies universality in the bulk for the m-point correlation 
function for the orthogonal polynomial ensemble (see Lubinsky [I7[ Section 

!])■ 

(b) The Theorem follows, for example, from [13, Theorem 1.1] combined with 
remark (e) there. The result was obtained previously, however, by both Simon 
[19] and Totik [22], using a modification of Lubinsky's approach from |16j . 

(c) The requirement of continuity can be relaxed to a Lebesgue point type con- 
dition, assuming boundedness and uniform positivity of w in a neighborhood 
of x. This is Theorem 1.2 of [18] . 

(d) Avila, Last and Simon [2], using a modification of the approach of [T8] , 
obtain universality for certain measures whose support is a positive Lebesgue 
measure Cantor set. 

(e) While it seems that local absolute continuity of \x is almost sufficient for 
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universality at x, it is certainly not necessary: Breuer jl] has recently shown 
there are purely singular measures such that universality holds uniformly for x 
in an interval (and a, b G R). Of course, K n (x, x) cannot be defined for purely 
singular measures, so the statement needs to be slightly modified (see jl] for 
details). 

In this paper we prove an analogue of Theorem 11.11 for averages fll.l[) . Here 
is our main result. 

Theorem 1.2. Let \x be a probabilty measure on R with finite moments. Let 
x G supp(/i) be such that: 

(i) There exists an interval, J, with x G J such that /x is absolutely continuous 
in J: dfi(t) = w(t)dt for t G J. Moreover, w G L°°(J, dt). 

(ii) w(x) > and x is a Lebesgue point of w, by which we mean that 

i r-x+r 

lim - / \w(t) — w(x)\dt — 0. 



(Hi) K n has a uniform universal limit at x. 

Under these assumptions, for any pairwise distinct aci 's and (3j 's, such that 
ai,..., a k el, Pi,...,0 k G C\R 



where 

(3 — a J ir(s — a)(s — (3) 

— oo 

As an immediate corollary of Theorems 11.11 and ll.2[ we obtain 

Corollary 1.3. Let \x be a probability measure on R with compact support 
that is regular in the sense of Stahl and Totik. Suppose x G supp(/x) has a 
neighborhood, J, such that fi is absolutely continuous in J: d(i(t) = w(t)dt for 
t G J, with w bounded in J. Assume further, that w is positive and continuous 
at x. Then for any pairwise distinct cti 's and (3j 's, such that a%, . . . , G R, 

A,...,AeC\M 



vi =i D n (x + ^7^) 

J \ K n (x,x) J 
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Remarks. 

(a) The function W(/3, a) can be rewritten as 

-«) 



Imp > 0, 



In. .»' < 0. 



/8-Q 

The formula above can be obtained from that of Theorem 11.21 by contour 
integration. 

(b) If d/j,(y) = const e~ v ^dy, and V(y) is an even polynomial, then Theorem 
11.21 reduces to the result obtained in Strahov and Fyodorov [21, Section 2.4]. 

(c) Note that the assumptions of Theorem II .21 do not include compact support 
of fi. The theorem essentially says that for any measure with finite moments, 
local absolute continuity combined with a uniform universal limit for K n at x 
imply uniform limits for averages (11 .ip at x. 

(d) As explained in remark (e) after Theorem II. 1\ there are purely singular 
measures for which K n has universal limits. As the absolute continuity of 
fi around x is essential to the proof of Theorem 11.21 we, unfortunately, have 
nothing to say about this case. 

The rest of the paper is devoted to the proof of Theorem 11.21 Namely, in 
Section [2] we present a useful algebraic formula for averages (II. ip (Theorem 
12. ip . This formula reduces the investigation of asymptotics of (II .ip to that of 
the Cauchy transform of the reproducing kernel, 

K n (t, a)djj,(t) 



Theorem 13.11 establishes universality for the Cauchy transform of the repro- 
ducing kernel. Theorem 11.21 (which is the main result of this work) is then a 
simple corollary of Theorem I2.1| and of Theorem 13.11 

2. A FORMULA FOR RATIOS OF CHARACTERISTIC POLYNOMIALS 

Theorem 2.1. Let 1 < k < n, and assume that a,\, . . . , . . . , fif. are 

pairwise distinct complex numbers. Moreover, assume that Im(fij) ^ for 
j = 1, . . . , k. Then we have 

where the two-point function W n ((3,a) is defined by 

p — a J t — p 

Proof. A formula of two-point function type for the average of ratios of charac- 
teristic polynomials was obtained by Baik, Deift and Strahov [21 Section III]. 
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Theorem 3.3 in |3j gives 
where 

~ /i n (/3)7r n _ 1 (o!) - h n -i(/3)ir n (a) 

W n [P,a) = . 

p — a 

In the formulae above 717(01) is the Zth monic orthogonal polynomial associated 
with /i, the function hi(f3) is the Cauchy transform of 7Ti(a), 

t hW) = 7 l r-j^m, 

2m J t — p 

and 

7W-i = -^ill-v 

Clearly, the formula for the average of ratios of characteristic polynomials 
can be rewritten as 

j-r D n (a>j) \ fc(*+D A(/3,a) A . nxr , a 

Ua^)/ =(_1) AfAK det(f " (ft ' a,)) «^' 

where 

7„_i /i n (^)p n _i(a) - /i n _i(/3)p n (a) 



7« /5 - « 

and where 

Pi(t)dfj.(t) 
t-p ' 

To obtain the formula for W n ((3, a) stated in the Theorem we note that 



7n-i h n (f3)p n - l (a) - h n -i{P)p n {a) 
In P - a 

7n-l 1 



p-a J t-P 



p n (t)dfi(t)\ ( f p n _i(t)d//(t)\ 
7— g ) Pn-i(a) - ( / : — 3 ) Pn{a) 

K n (t, a)dfi(t) 



P — a \\J t — P J \J t — P ) 

t — a 



1 f K n (t,a)dp,(t) 



P-a J t-P 

where in the second equality we have used formula (jl.2p . and in the third 
equality we have used the reproducing property of K n . □ 
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3. Universality for the Cauchy transform of the reproducing 

KERNEL 

Theorem 3.1. Let x G supp(/x) be such that conditions (i) — (Hi) of Theorem 
are satisfied. Then for any a6l, /3 G C \ R 



lim — 

n ^°° K n (x, x) 



\ K n (x,x)' 
t — X 



dn(t) 



where 



S(a>, s) 



K n (x,x) 



sin TV (a 



at, s. 



ft 



ds, 



Tv{a — s) 

Remark. In case of the Chebyshev weight for example, 

1 



dfi(y) = w(y)dy, w(y) 



, ye (-1,1), 



Theorem 13.11 can be checked by direct computations. These computations are 
similar to those of [T7J Section 2], for the reproducing kernel associated with 
the Chebyshev weight. 

Theorem 13.11 follows from the following 

Lemma 3.2. Fix a G K and G C\R. Under the conditions of Theorem \l.S\ 

for any M >2\(3\, 



lim sup 



(3.1) 



K n [x + 



K n (x,x) ' 



< 



■K-ti (x, X 

S(a, s 



t — x — 



dfi(t) — 



K n {x,x) 



S(ot, s) 
s-P 



ds 



[-M,M] 



s-P 



■ds 



+ 



\w\\j 
w(0)M 



where \\w\\j is the essential supremum of w on J. 

Proof of Theorem \3.1\ assuming Lemma \3.B. Fix a e I and ft G C \ R. As 
M — > oo the righthand side of inequality (13. ip goes to zero. This is because 



sin (tv is 



a)) 



s — ft tv(s — a)(s — ft) 

is integrable so its tail goes to zero. Since M is arbitrary, we get the limit 
relation in Theorem 13.11 □ 

Proof of Lemma \3. M By shifting the measure, we can assume that x — 0, 
which we henceforth do for ease of notation. Fix a G R and ft G C \ R, and fix 
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M > 21/31. Write 



K n (0,0) 
1 



n V^n(0,0) 



d/i(t) 



K n (0,0) 



8-0 



ds 



< 



K n (0 } 0)J t- 



h " ' £„(0,0)'* 



P 



dn(t) - 



M 



M 



§(«, s) 

~7^p 



ds 



K n (0,0) 

To prove the Lemma, it is enough to show that 



-M,M] 



S(a, s) 
s-P 



ds 



lim sup 



(3.2) 



1 



K n (0,0) 



P 



d/i(t) - 



K n (0,0) 



M 



M 



s-P 



ds 



< 8i 



\w\\j 
w(0)M' 



Let I n = [- 



W^ andwrite 



i /H 

'0,0) J t- 



V^n(0,0)' 



K n (0,0) 
1 



P 



dji{t) 



K n (0,0) 



< 



+ 



K n {0 } 0) Ji n t- w 



(0,0) 



dn(t) - J 



M S(a, s) J 
W ds 

—M S-P 



M 



8-fi 



■ds 



1 



K n (0,0) 



* - 



dfi(t) 



K n (0,Q) 

Consider first the first term in the righthand side of the inequality above, 



n \K n (o,oy 



K n (0,0)Ji n t- 







d [/,(£) — 



K n (0,0) 



M 



M 



h{a, s) 
s-P 



ds 



Since is not a pure point of [/,, it follows that K n (0,0) — > oo. Thus, for 
sufficiently large n, I n C J. We can therefore write 



K n (0,0) Ji n 



K, 



( — 

11 \K n {0,0)> 



t- 



P 



K n (0,0) J i, 



dfi(t) - J 



M S(a,s) 



K n (0,0) 



M 



s-P 



ds 



t - 



. P 

K n (0,0) 



/M 
-M 



~s~P 



ds 
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and by transformation of variables t 



K n (0,0) 



we get 



K. 



K„(0,0) 



K n (0,0)Ji n t 



,t)w(t)dt r M 

M S P 



K n (0,0) 



i 



ds 



M 



M 



K n (0,0) ' K n (0,0) 



w 



K„(0,0) 



8-0 



K n (0,0)w(0) 



— S(a, s) ds 



The righthand side of the equality above goes to zero as n — > oo. This follows 
since K n has a unifrom universal limit at 0, since w is essentially bounded on 
J and since is a Lebesgue point of w (note that \s — /3\ > | Im(/3)| > 0). 
It remains to show that the inequality 



(3.3) limsup 



K n {o,o) 



■t 



K n (0,0)Mi n t- 



d/i(t) 



< 



K n (0,0) 

holds. First, by the Cauchy-Schwarz inequality we obtain 



\w\\j 
w(0)M' 



K n (0,0) 



dfi(t) 



< 



K n (0,0) 



( 



V 



K n (0,0) 

dji{t) 
(t-S—\ 

V K n (0,0)J 



1/2 



1/2 



a 



I 



^„(0,0) 



,* dn{t) 



Second, the above inequality, and the reproducing property of the kernel imply 



— I 

(0,0) Jr 



K n (0,0) Jm\i n t 



K n (0,0) 



< 



K ( ^-2 ^ 

n \K n (0,0)> K n (0,0) J 



1/2 



w(p)K n (p,0) 
Note that for sufficiently large n, 



1/2 



dfi(t) 



(t--e—\ 



a a 



tf„(0,0)' K n (0,0) 



1/2 



w(0)K n (0,0) 



< 
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by assumption (Hi) in the statement of Theorem 11.21 Thus we get, for suffi- 
ciently large n, 



(3.4) 



1 



K n (0,0) 



K„ 



K n (0,0) 



■t 



t 



8 



K n (0,0) 



< 



w(o)K n (o,oy/ 2 

We are left with estimating 



dfi(t) 



djj,{t) 



1/2 



t — 13 

K n (0,0) ) J 



djj,{t) 



t - 



8 



2 ' 



First note that, since M > 2| 



K n (0,0) 

we get for every t with \t\ > 



M 



tfn(0,0)' 







Therefore, 



K n (0,0) 
dfi(t) 



t 

> -. 

~ 2 



< 



t 



Now we split the integral in the righthand side of the inequality above. Let 
Hn = [~ x n( o o)i/3 ' K ni o 0)1/3 ]' and a S ain note tliat for sufficiently large n, H n C 
J. Clearly, I n C H n , so we can write 

R\I n = (R\H n )U(H n \I n ). 

We split the integral accordingly, 



4d/z(t) 



t 2 



< 



Ad/i(t) 
4:d/i(t) 
Adjx(t) 



4d(i,(t) 



H n \I n 



t 2 



+ i 

H n \I n 



Aw(t)dt 



t 2 



dt 



H n \In P 



where we have used the fact that H n C J to write dfi(t) = w(t)dt there. The 
proof of Lemma [3.21 is almost complete. For the first integral on the left, note 



that for t $ H n ,t 2 > 



M- 



. Therefore, we have (taking into account that 
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H is a probability measure 



t 2 



< 4 



M 2 



2/3 



djj,(t) = 4 



ifn(0 ; 0) 2/3 

M 2 



For the second integral, integration of 4 gives (recall the definition of I n ) 



A\\w\ 



H n \I n 



dt 



< A\\w\ 



dt 



\ w \\j- 



M ' 



Plugging these estimates into (13.41) . we see 



K n (0,0) 



if„(o,o)' 



t - 



8 



■dfj,(t) 



(3.5) 



< 



w(0)X n (0,0)V2 



K n (0,0) 

K n (0,0) 2 ^ 
M 2 



+ \\w\\j 



K n (0,0) 



1/2 



M 



1 II II 1 



1/2 



which immediately shows that 



lim sup 



K n (0,0) 



K„(0,0)' 



t - 



dfi(t) 



K„(0,0) 




w(0)M' 



□ 
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